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Extending previous results [Phys. Rev.D 56, 6351 (1997)], we estimate the back reaction of 

cosmological gravitons in the expansion dynamics, during the matter age. Tensor perturbations 

with scales larger than the Hubble radius are created due to the inequivalence of vacua at different 

times. During noninflationary phases these perturbations become effective gravitational waves as 

they enter the Hubble radius, adding new contributions to the energy density of the subhorizon 

waves. During the radiation epoch the creation of these effective gravitons may lead to a departure 

from the standard behaviour a{t) oc t^'^, with possible consequences to several cosmic processes, 

such as primordial nucleosynthesis. We examine the implications of this phenomenom during the 

. . . matter dominated era, assuming an initial inflationary period. The dynamical equation obeyed by 

CTS I the scale factor is derived and numerically solved for different values of the relevant parameters 

O^ . involved. 

0^ 






oo 
en 



^ 



PACS number(s): 04.30.Nk, 04.30. Db, 04.62. +v, 98.80. Cq 



I. INTRODUCTION 



In a recent paper |y] , we have studied the back reaction of cosmological gravitons on the cosmic dynamics during the 

radiation epoch. The analysis was based on the foUowing reasoning: the inequivalence of vacuum states at different 

^ moments of time leads to the production of tensor perturbations in scales larger than the Hubble length [1 - 8] . During 

noninflationary periods of expansion, these very long tensor perturbations (VLTP's) become effective gravitational 

1^ _ waves (EGW's) as they enter the Hubble radius, thus adding new contributions to the energy density associated with 

fvq . the subhorizon waves pg [1-5], [|j. Such a phenomenum can be studied as a process of creation of effective gravitons 

f^ ' by using the macroscopic formalism to matter creation based on the thermodynamics of open systems |10| , |ll| . In 

On • order to deal with the process, a creation pressure term is introduced in the continuity equation obeyed hy pg. A 

On ' dynamical equation for the scale factor a{t) is then derived which takes into account the effective gravitons' back 

r-| \ reaction. Full details of this approach can be found in Ref. 0] where the equation for a(t) was numerically solved for 

^ i' a model in which the universe evolves from an arbitrary initial phase to a radiation dominated era. It was found that, 

' if the barotropic index 7 of the equation of state in the first epoch is close to 2/3, the back reaction of the effective 

^ gravitational waves makes a{t) deviate appreciably from the standard behavior a{t) ex t^". 

"^ ■ In the present paper we extend our analysis to the matter dominated period of the cosmic expansion. In particular, 

C^ , we aim to check a conjecture stated by Sahni [ [L2[ , according to which the process mentioned above would ultimately 
ILj^ ' lead the universe to expand linearly with time (a(i) oc t), since the tensor perturbations would enter the Hubble 
radius during noninflationary periods of expansion (increasing pg) and leave the Hubble radius during inflationary 
periods (decreasing pg). In Ref. [Q we have shown that the scale factor would turn from a{t) oc t^/"^ to a{t) oc t only 
if 7 ^ 2/3. Several models with the universe evolving linearly with time during its early phases have been proposed, 
most of them related to string motivated cosmologies |l^,|lj] . An universe evolving hnearly with time during most of 
its history would lead to important observational consequences, such as those related to the age of the universe, to 
the luminosity distance - redshift and the angular diameter distance - redshift relations, and to the galaxy number 
count as a function of the redshift [fsUlq] (see also |17[]). 

We should remark that, in order to the mechanism of graviton creation take place, no 'exotic' physics has to be 



assumed. It only requires the validity of quantum mechanics and of general relativity 1 18|. Therefore, the study of the 
possible influence of the cosmological gravitational waves in the expansion dynamics has an importance on its own, a 
fact that has not been fully appreciated in the literature. 

The paper is organized as follows: In Section||we present the basic equations related to the spectrum of cosmological 



gravitational waves. In Section IH we derive the dynamical equations that govern the expansion dynamics during the 
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matter age, taking into account the gravitons' back reaction. In Section IV we show the numerical solutions of these 
equations and present our conclusions. 

The system of units used is such that h = c ^ ks ~ I- 

II. BASIC EQUATIONS 

We will consider a homogeneous and isotropic universe, for which the background line element takes the Friedmann- 
Robertson- Walker (FRW) form 

ds^ = dt^ - a^{t)df = a^{f]){dri^ ~ dP) , (1) 

where t and rj are, respectively, the cosmic and conformal times, related by 

dt = adrj . (2) 

We will restrict ourselves to the spatially flat case. If the pressure p and energy density p of the cosmic fluid are 
related by the equation of state 

P=il-l)p, (3) 

then the scale factor is H 



a{t) = ao 



'-^(t-to) 



2/(37) 



(4) 



where Hq = H(to) and H = d/a. 

We will further assume that the universe evolves from an initial arbitrary era (7 = 7) to a radiation dominated 
phase (7 — 4/3) and then to a matter dominated period (7 = 1). 

As the inequivalence of vacua appear at different instants of time, we will focus on the time-dependent amplitude 
of the tensor perturbations fj,{k,r]), which obeys ||l|,||,|| 

^,'\k,v)+le-^W{k,r^)^0. (5) 

In the above equation the primes indicate derivatives with respect to the conformal time and k is the comoving wave 
number, related to the physical wavelength A and frequency uj by 

27ra 
k — — — = uja . (6) 

A 

The properly normalized solutions of Eq. (|5|), corresponding to a scale factor given by Eq. (j^, and representing 
adiabatic vacuum states |0,||, can be found in Refs. [^|J^. 

The spectrum of the EGW's can be described by the quantity Pg{iLj), defined in such a way that Pg{uj) dcu represents 
the energy per unit volume between the frequencies w and oj + duj [^-Q . Note that this quantity can be defined only 
for the EGW's 1^,^, i.e., for those perturbations such that 

\<\h = H-K (7) 

In the units used in this paper and assuming an initial vacuum state, Pg{uj) is given by |0-0| 

Pgicj) = ^ (7V(^)) . (8) 

In the above equation {N{uj)) is the expectation number of the gravitons with frequency uj which is found to be [|l],Hj] 

{Niu;)) = |/3(fc)p , (9) 

where /3(fc) and a{k) (see below) are the Bogoliubov coeficients relating the creation and annihilation operators 
(which define the particle states) of different epochs. We will denote by ai, (3i the coeficients associated with the first 



transition from an arbitrary phase to the radiation dominated era, and by a2, P2 those related with the transition 
from the radiation to the matter dominated epoch. These coefficents have been obtained in references [g|Q] (see also 
0). The coefficient (3 is then given by ||] 

/3 = /32(fc)ai(fc)+a^(/c)/3i(fc), (10) 

and the asterisk indicates the complex conjugate of a quantity. 

The total energy density associated with the tensor perturbations is obtained by integrating Pg{Ld), 



Pg{t)= PgiL^)du, (11) 

where 

Urmnit) = 27riJ(i) . (12) 

This infrared cutoff is obviously related to the condition (H) and, since it is time dependent, it is the origin of the 
process of creation of EGW's. In a FRW scenario, it will be responsible for the fact that pg does not decay with a~'^, 
as one could expect for a massless particle such as the graviton. 

The integral in Eq. (|l^) is usually written with a finite upper limit. This is done because the above method for 
obtaining the spectrum, based on the sudden transition approximation between cosmic eras, does not give accurate 
results for high frequency waves. In fact, it is reliable only for those waves whose periods are much greater than 
the transitions time scales |4|jl9|]. However, due to the adiabatic theorem ||], the number of created particles should 
decrease exponentially for the high frequency modes. Hence, an ultraviolet cutoff is imposed and the error produced 
by doing so is supposed to be small. Our present results will not be affected by this approximation since we are 
obviously interested only on the very long modes. 

III. DYNAMICAL EQUATIONS DURING THE MATTER ERA 

Following Rcf. |l|], we will suppose that, prior to the time ti the dominant material content of the Universe has an 
equation of state of the form (H) . This stage can be cither an inflationary or a noninflationary one. At ti a transition 
occurs, so that the new barotropic parameter is 4/3. Due to the transition, tensor perturbations are created ||l|,0,D. 
Some of these have wavelengths less than H~^{ti) (EGW's), but most are created with scales larger than H^ {ti) 
(VLTP's). A similar phenomenon! will occurr at a time ^2 when the universe becomes matter dominated (7 = 1). 
We will further assume that during this matter dominated period and until a time, say, tg, the energy density pg 
associated with the EGW's is negligible compared with the energy density of matter p,„ . There is no restriction over 
the size of the interval ig — ^2, that can be taken to be arbitrarily small. However, after tg^ the ongoing transformation 
of VLTP's into EGW's makes the continuous increase of pg start perturbing the dynamics. Thereafter, the Universe 
can be supposed to be filled by two fluids: dust and the one composed by the effective gravitons. The total energy 
density is then pm + Pg- As gravitons require extremely high energies to interact with matter P,pO[, the two fluids 
can be safely considered to be noninteracting. The effective gravitons behave as a perfect fluid with equation of state 

p. = f . (13) 

Nevertheless, it is possible to associate a creation pressure term Hg to the creation process of these gravitons Q. The 
field equations are then written as 

(14) 



• 2 

8lTG{pm + Pg) = 3— , 
0/ 




Pm+Pg+Ilg)^-2-- 


d2 


rig + ZHUg = ^g , 





8^G(p™+Pg+ng) = -2--— , (15) 

(16) 



where the dots indicate derivatives with respect to the cosmic time, Pm = is the pressure of the dust fluid and Ug 
and \l/g are the number density and creation rate of the effective gravitons, respectively. This last equation is the 
novel aspect introduced by the phenomenological formalism to particle (in this case, gravitons) creation |l^Jll|] . The 
conservation of the energy-momentum tensor leads to 



where p is the total energy density [p = p,,„ + pg) and p is the total pressure [p = Pm + Pg, Pm = 0). 
As the two fluids are nonintcracting, this conservation equation can be split into 



p',n + iHprn = 



and 



Pg + iH\ -Pg + n 



From Eqs. (|l|), (|l|), ® and (|l|) we obtain § 



a 1 a^ 

-+o- 
a la 



7+0^2 =4^G(p, + ^^ 



iH I ' 



(17) 



(18) 



(19) 



(20) 



a2 SvrG 



(Prn + Pg) ■ 



(21) 



It is important to mention that, in the method described in |0,0,B, the Bogoliubov coefflcients that appear in the 
calculus of (iV(u;)) are evaluated at the transition times ti and t2. Hence, Pg{t) is univocally determined in terms of 
a(t) and d(i), that is Q] 



Pgit) = Pg{a{t), a{t)) . 



(22) 

Therefore, the system of coupled equations (|l8|), (|2l| ) and (^ allows us to obtain a(t) taking into account the back 
reaction induced by the transformation of very long tensor perturbations into effective gravitational waves. 
Let us define 



A{t) 



a(i) 

02 



(23) 



and 



02 = 0(^2) , 



(24) 



oi = a{ti) , 



(25) 



Hi = H{h) , 



(26) 



02 

Ol 



(27) 



Then Eqs. (ph-(0), (ttq), ([2l|), and the expressions derived for the Bogoliubov coeflcients in P,01, lead to a lengthy 
set of equations that enable us to evaluate a{t) during the matter age. These equations take the form 



^ - o^ A 



\mpij A 



(28) 



where mpi is the Planck mass. The values of the constant k and of the function G{t) depends whether 7 = (the 
initial era is a de Sitter one), or 7 7^ (we will restrict our analysis to inflationary models, so that 7 < 2/3). 
For 7 = 0: 



Git) = 1 + 





2 






'^=3^' 




1 ^ 

16a* ^ 


Hi 

T = 


ybi 



+ C0 , 



b2a^A 



(29) 
(30) 

(31) 



Co = 167r 



1 1 



bj bta^J ' 



(32) 



jit)=Ni{(T, T)s[n9i+ N2{(T, t) cos 9i + Ns^a) cos 02 + Ni{a) sin 02 
+ N5ia)J{t), 



(33) 



9i(i) =47ra(cr-l) 



A 
^ 



(34) 



47r(cr - 1) 

62 CT 



(35) 



r=^ cos 6* 

"'ei(t) ^ 



(36) 



Ni{a,T) 



A-K 



(3cr+l) 3 62(llcr3 + llcr2- 7(7 + 1) 



(37) 



^2(a,T)EEE- 



(110-2 +6(7-1) 



^2 "2 _4 



(38) 



^3(^) ^ 2 



■2 (11(7^ +6(7-1) 

8^ a^ 



(39) 



^^(-) - fe 



fo2(llo-3 + 110-2 -7(7 + 1) (3cr + l) 



(40) 



iV5(^) = - 



- ^2(0-2) (^^^3 + 11^2 _ 7^ ^1) 



(41) 



On the above equations the free parameters 61 and 62 are related to the transitions time scales Ati and At2 that 
are assumed to take place at ti and ^2, respectively: 



Ail 



bi _ bi 
H{ti) Hi ' 



(42) 



Ai2 = 



62 _ 62 

mr)^H~2 



(43) 



Thus, the larger bi, the slower is the transition in comparison with the Hubble time at U, H~ . It is usual to take 
&i ~ 62 ~ 1, but, as for the first transition, there is no compelling reason to assume that it necessarily ocurred in a 
time scale of the order of HT in every model. For example, this transition may occurr as fast as 10^^ H^^ in the 
'new' inflationary scenarios ||]. 

When the first stage is not a de Sitter one, that is, for 7 < 2/3, 7^0, the constant k and the function G{t) become 



4 

K = — 

3 



(44) 



where 



G{t) = -^Giit) 



TTCT"' 



\2m + 



Y^Co, 



3 (2-7) 
2 (37 - 2) 



(45) 



(46) 



Co = \(m-l 



m + - 



yi 



B„,{y)dy + D{yi) ~ D{y2) 



V2 



(47) 



J/i 



y2 



7r|2m + 1| 
bl 



7r|2m+l| 
620- 



Diy) = fiim,,y)Bmiy) + f2im,y)Bm+iiy) + f3im,y)Crn{y) 



(48) 
(49) 
(50) 



fi{m,y) = 



^'+U'"+^)('"+^^^' 



--\m-- \m + -]y 



(51) 



/2(m,y) = - 



— \ rn ^ — 
2 V 2 



^^-^ly 



(52) 



h{m,y) = "2 ("^+ 2 



v' + \[^-\] v' 



(53) 



B^{y)^Jl{y) + Y^{y), 

Cm{y) = Jm{y)Jrn+l{y) + Y.m{y)Y.^+i{y) . 

Jm and Y„i are, respectively, the Bessel functions of the first and second kinds of order m, and 

1 



Gi (i) = -- In T + P(m, (7, 62) 



62 



53(-2m-3) 
1 



(^ 



-2m-3 



1 



5127r2(l-2m) 



/ l-2,n _ ^^ 



1287r 



/W 



(54) 
(55) 



(56) 



P(m, cr, 62) = 



2 6- 



2Cr 



TT 2to+ 1 



-2m 



r2(-rn)|2m+ll 



(57) 



/(i) =Fi(m,cr)r-2™sin6'i 

^2 _l-2m 



47r(l - 2m) 

Pi (m, ct) sin 6*2 + P4 (m, cr)/(i) 



^P2K^)r-2'"-^- 
02 

cos 01 + Pz{m, a, 62) cos 6*2 



smi 






Pi{m,a) 



P2Ka). ("-^) 



cr2|2m+l| 



3cr-2m(l + 2cr) 
(-2m)(l-2m)CT ' 

3(7 - 2m(l + 2cr) 2cr(CT + l) 
(--2m)(l-2m) ^ (cr - 1) 



P3(m, cr, 62) 



47r(l - 2to) 62 



47r 
-—P2{m,a) 



(58) 
(59) 

(60) 
(61) 
(62) 



P4(m, 0-) 



and r is the gamma fmiction. 



1 



4(a - 1) 



a\2rn+ 1 
2(T2|2m + l 
*" (a -1)2 



|2m+l| 



3cr - 2m(l + 2cr) 2cr(cr + l) 
(-2m)(l-2m) ^ (ct - 1) 



(63) 



IV. NUMERICAL RESULTS AND CONCLUSIONS 

We are interested in evaluating a{t) after ^2, hence the initial conditions are 

Aih) = 1 
and 

A{h) = ^ . 

Note that, besides 61 and 62, there are three free parameters: 

hi = Hi/mpi , 
(7, whcih we relate to the duration of the radiation era i2 ~ ^1 through 

C72 = l + 2i7l(i2-il), 



(64) 



(65) 



(66) 



(67) 



and m (or, equivalently, 7). 

We have integrated Eq. ( |2q ) numerically for different values of the parameters &i, 62, ^1, 7, and ^2 — ^i- We have 
taken as characteristic values 61 = 1, 62 = 1, ^1 = 10^^° and < 7 < 0.6 . In Figure 1 we show the solutions for five 
values of ^2 — ti, namely 3 x 10^ yr (curve a), 10^ yr (curve b), 3 x 10^ yr (curve c), 10^ yr (curve d), and 3 x lO^'yr 
(curve e) . The time is measured in units of the time of the beginning of the matter era ^2 ■ The important thing to 
notice is that the solutions are insensitive to all parameters, except ^2 ^ ti- We see that different values of ^2 — ii 



change the behaviour of A(t) near ^2- All the curves go assymptotically as t^^^, which is represented by the dotted 
curve. However, the net effect of decreasing the duration of the radiation era is to change the present value of a{t)/a2. 

We are led to conclude that the transformation of very long tensor perturbations into effective gravitational waves 
does not change the expansion dynamics during the matter epoch. This is in contrast to what happens during the 
radiation era, as it has been shown in Ref. M. Therefore, the Sahni's conjecture |12] will hold only under very 
restrictive conditions. 

We must remark, however, that, in the model analysed above, we have assumed a standard evolution from a radiation 
phase with a{t) ex i^/^ to a dust era with a{t) oc i^/'^. We have not taken into account the modifications in the dynamics 
generated by the gravitons' back reaction during the radiation dominated period M. A more realistic model would 
consider these modifications and then a transition to a matter dominated era. The numerical work becomes much 
more involved in this case. Moreover, primordial nucleosynthesis and the data related to the anisotropy of the cosmic 
microwave background will place severe constraints on the relevant parameters. These questions are presently under 
investigation. 
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FIG. 1. A logarithmic plot of the quantity A = a/a2 as a function of time [solution of Eq. (28)] for bi — 1, b2 = 1, hi — 10"^", 
7 = 0.5, and five values of i2 — ti: 3 x 10® yr (curve a), 10® yr (curve b), 3 x 10^ yr (curve c), 10^ yr (curve d), and 3 x 10'* yr 
(curve e). The dotted curve corresponds to a{t) oc t^''' and the log is to base 10. The solutions are insensitive to all parameters, 
except t2 —t\. 



